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1 Introduction 

Recently it was reported, that the polarity of the T- 
wave in magnetocardiograms (MCG) recorded 
during pharmacologically induced stress changes 
with respect to rest in a subset of the MCG channels. 
This often leads to opposite R-T polarity. Contem¬ 
porarily it is reported, that no comparable changes 
can be observed in ECG-recordings [1]. This finding 
could be the first experimental evidence for a dif¬ 
ferent information content of magneto- and electro¬ 
cardiogram. 

In the very first report on cardiomagnetism [2] it 
was suspected that closed loop currents within the 
heart might set up detectable magnetic fields while 
producing no detectable electric potentials. It was 
theoretically shown in successive studies [3,4] that 
in structures of anisotropic conductivity the non¬ 
vanishing curl of the current density will result in 
the generation of closed current paths (vortex 
source) giving rise to electrically silent magnetic 
fields. However, for a complex three-dimensional 
fibrous structure like the ventricular muscle, the 
splitting of the current density field in a vortex part 
and a source-sink part is a computationally deman¬ 
ding task. Thus, up to now no theoretical estimation 
exists at which order of magnitude electrically silent 
magnetic fields can be expected. 

This study experimentally and theoretically in¬ 
vestigates magnetocardiographic T-wave inversion. 
Measured data has been obtained in two healthy 
male subjects. This data can be compared to 
simulated field patterns obtained for a simplified 
model of the left ventricle embedded in a bounded 
torso model of homogeneous conductivity. An 
estimation of the magnetic field components 
produced by closed loop currents is obtained by a 
dipolar approximation of the source (current dipole 
and magnetic dipole). 


2 Methods 

2.1 Measured data 

Two male subjects, 33 and 35 years old and both of 
them cardiologists, volunteered for stress MCG 

recordings. Measurements were carried out at the 
Biomagnetic Center, Jena, Germany. Similarly as in 
[1] stress was induced pharmacologically in order to 
avoid involuntary movements of the subjects 

relative to the sensors. A stress test protocol was 
applied which is standard for stress echocardio¬ 
graphy investigations. First, a continuous injection 
of Dobutamine was given in steps of 5, 10, 15, 20, 

30 and 40pg/kg/min. Each step lasted 3 minutes. 
Finally, a bolus of lmg Atropine was given in one 
subject and 0.5mg Atropine in the other in order to 
achieve a maximal heart rate of about 160 beats per 
minute. 

Magnetocardiographic measurements were carried 
out in an actively and passively shielded room, 
which ensured a low level noise floor for both low 
and high frequencies [5]. A dewar consisting from 

31 axial SQUID-gradiometers (Philips, Hamburg) 
was positioned about 2cm from the left anterior 
chest wall. Care was taken, that about the half of the 
signals had positive amplitudes and the other half 
negative. 

The amplifier system used for the acquisition of 
electric potentials (Philips, Hamburg) allows a 
maximum of 32 channels. Thus, body surface 
potential (BSP) mapping was carried out on the 
anterior chest wall only. A mapping scheme pro¬ 
posed in [6] was adopted (‘Lux-anterior’). The zero- 
potential was defined by a Wilson-terminal. 

2.2 Numeric field calculation 

A bidomain model based BEM-FEM coupling 
technique is applied in order to properly consider 
myocardial tissue properties for the calculation of 
the electric potentials and the extracorporal 
magnetic field [8]. Here, finite elements are used for 
modeling the anisotropic myocardium. In the 



surrounding homogeneous and isotropic volume 
conductor the boundary element method (BEM) is 
adopted. The zero potential is defined in such a way 
that the sum of the potentials at three nodes located 
in the lower left abdominal region and at the right 
and the left shoulder gives zero (Wilson terminal). 
Applying BEM-FEM coupling the electric potential 
(p can be calculated in each point within the ventri¬ 
cular volume Q v and in each point at the surfaces of 
torso and ventricle F T and T v . The extracorporal 
magnetic field can be calculated applying Biot- 
Savart’s law. For the isotropic torso region 
integration can be reduced to a surface integral [8,9]. 
The current dipole moment p associated with the 
distributed source is given by: 

p = - Jcr. gradV m d£l, (1) 

where V m denotes the transmembrane potential 
(source term) and <7/ is the effective intracellular 
conductivity [ 8 ]. 

The magnetic dipole moment m associated with the 
distributed source can be calculated by [9]: 
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where r f denotes the vector from the source point to 
the location of the dipole and n is the surface normal 
pointing out of the ventricular mass. The bulk 
conductivity of the ventricle is named k b and k t is 
the conductivity of the torso region. Here, the choice 
of the dipole location enters the result. Also the 
electric potential (p enters the calculation in addition 
to the independent source term V m . The location of 
the magnetic dipole was assumed to be at the center 
of gravity of the ventricle. This choice ensures that a 
homogeneous current density field in the ventricles 
will produce almost vanishing magnetic moments. 

2.3 Model description 

An ellipsoidal model of the left ventricle was placed 
in realistic position in a realistically shaped torso 
model. The ventricular model was adopted from [7] 
where it is described in detail. In short, the 
ventricular volume L2 V can be described by the 
curvilinear coordinates 77 , 7 ?, and r: 

x v = a (r)cos $ cos 77 \7i <$<2n 

y v =a(r)cos$sm rj 0< 77 < 2zr , ( 3 ) 

z v = c{r ) sin $ 0 < r < 1 


with 

a(r) = a x + r(a 2 - a x ), c(r) = c x + r(c 2 - c x ) . (4 ) 

Here x V9 y v and z v form a local Cartesian coordinate 
system. The parameters a 7 , a 2 , Cj and c 2 define the 
geometric dimensions of the ventricle. For the 
endsystolic geometry at physical rest they have been 
assigned to: aj = 1.5cm, a 2 = 3.4 cm, c 7 = 4.7 cm and 
c 2 = 6.0cm. Stress was modeled by an increased con¬ 
traction, keeping the muscle volume almost 
constant. This was done by assigning a 2 = 0.7cm, a 2 
= 3.1cm, Cj = 3.4cm and c 2 = 4.7cm to the para¬ 
meters. Cardiac fiber orientation is identical as sug¬ 
gested in [7] rotating in counterclockwise direction 
from the epi- to the endocardium. When positioning 
the ventricular model within the torso (transforma¬ 
tion to a global Cartesian coordinate system x 9 y and 
z) care was taken, that the location of the heart base 
remains almost the same (position of the heart is 
fixed by the great vessels attached to the atria). 
During repolarization the M-cells yield a prominent 
contribution to the body-surface-potential (BSP) [ 8 ]. 
For the time instant of the T-wave maximum the 
transmembrane potential has the lowest value at the 
epicardium and the highest value in the mid¬ 
myocardium. The transmembrane potential at the 
endocardium is in-between. The transmural 
variation of the source term is modeled by a radial 
dependency of second order. Additionally, a small 
azimuthal variation is included into the model 
considering transmembrane potential differences in 
the free ventricular wall and the septal wall. 

F>,^) = F 0 +AF r (f + fr-|r 2 ) 

( 5 ) 

AV 

+ y [l - sin $ sin(77 + )] 

Here the following values are assigned to the para¬ 
meters: V 0 = -90mV, AV r = 60mV, AV r = 4mV and r/o 
= -k!9. For the anisotropic model the following 
values are assumed for the bidomain conductivities: 
longitudinal bulk conductivity 0.60Sm _1 , transversal 
bulk conductivity 0.158nf 1 . longitudinal effective 
intracellular conductivity 0.30Sm _1 and transversal 
effective intracellular conductivity 0.03Sm''. For the 
isotropic model the bulk conductivity is 0.458irf 
and effective intracellular conductivity is O.OSSm ' 1 
[ 8 ], The conductivity of the torso domain is assumed 
to be 0.20Sm _1 . The left ventricular mass is modeled 
by 6240 linear thetrahedral finite elements (1605 
nodes). 1400 triangles on the boundary of the ven¬ 
tricle are used for BEM-FEM coupling. The torso 
surface is modeled by 1090 triangles (547 nodes). 
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Results 


Simulated data 


3.1 Measured data 


In both subjects T-wave inversion was observed in a 
subset of the MCG-channels, while no comparable 
changes were observed in the ECG. Fig. 1 shows 
sample MCGs. MCG T-wave inversion was already 
established at a heart rate of about 90 to 100 beats 
per minute (bpm). No further changes in T-wave 
polarity were observed with increasing rate (up to 
160 bpm). Thus, as Dobutamine is known to have a 
strong positive inotropic effect one can suspect that 
T-wave inversion is associated with changes of 
cardiac geometry due to increased contraction. 

Fig. 2 shows the magnetic field maps (MFMs) at 50 
and 90 bpm for the time instant of T-wave maxi¬ 
mum in Wilson-lead V 5 . Sensors located in-between 
the isofield-zero-lines at rest and during stress dis¬ 
played T-wave inversion. Movement of the subject 
relative to the sensor was less than 3mm as could be 
verified by 5 location coils placed on the anterior 
chest wall. 
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Figure 1: MCG (raw data) measured by the sensor 
at the center of the gradiometer array at a heart 
rate of 50 bpm (left panel) and at 90bpm (right 
panel). Note that T-wave inversion leads to opposite 
R-Tpolarity during stress. 



Figure 2: MFM during T-wave maximum at a heart 
rate of 50 bpm (leftpanel) and 90 bpm (right panel). 
Isofield-lines are shown in steps of 2pT. Dashed 
lines indicate negative values. The bold lines define 
zero. Note the remarkable displacement of the zero¬ 
line. This was emphasized by showing the zero-line 
at 50 bpm also as a dotted line on the right (see also 
text). 


3.2 

The BSP maps and MFMs computed for the models 
of anisotropic conductivity are shown in Fig. 3. The 
MFM-zero-lines displayed a remarkable displace¬ 
ment comparing the maps computed for rest and 
stress. The isopotential-zero-lines showed no com¬ 
parable dislocation. However, the simulated MFM- 
displacement had almost opposite direction with 
respect to the measured. 

For the model of isotropic conductivity no zero-line 
displacements were observed for both the BSP map 
and the MFM. 




rest stress 



Figure 3: BSP maps (upper row) and MFMs (lower 
row) computed assuming anisotropic conductivity. 
The torso is depicted in anterior view. Isopotential¬ 
lines are shown in steps of lOOpV. Dashed lines 
mark negative potentials. The bold full line marks 
the isopotential-zero-line. Isofield lines are drawn in 
steps of 2pT. Coding is the same as used for 
isopotential-lines. The zero-lines computed at rest 
are also shown as dotted lines in the maps computed 
for stress making displacements visible. 

Current and magnetic dipole moments were com¬ 
puted in order to provide a coarse splitting of the 
flux and vortex part (see Table 1). As they are of 
different physical dimension, they can not be 
compared directly. Thus, the magnetic fields pro¬ 
duced were computed for comparison. It turned out 
that for the model of anisotropic conductivity the 

















vortex components of the magnetic field are about 3 
(rest) to 5 (stress) times smaller with respect to the 
flux components. For the isotropic model the vortex 
components were always more than hundred times 
smaller than flux components. 

Table 1: Current and magnetic dipole moments 
computed for the anisotropic (a) and the isotropic 
(i) models, x is defined from right to left , y from 
anterior to posterior and z from caudal to cranial. 
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nAm 2 
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a, rest 

2.9 

-2.1 

-1.6 

20.3 

-5.5 

-20.5 

a, stress 

2.3 

-1.9 

-1.1 

9.0 

-3.0 

-8.7 

i, rest 

7.3 

-2.4 

-6.5 

0.9 

-0.6 

-1.6 

i, stress 

5.7 

-2.0 

-5.1 

0.9 

-0.4 

-1.4 


4 Discussion 

One goal of this study was to investigate the role of 
closed loop currents in the genesis of T-wave inver¬ 
sion in MCG stress recordings. The simulated 
magnetic dipole moments are in the same order of 
magnitude as those estimated inversely in [1]. 
Simulations showed that vortex components of the 
magnetic field should be detectable in principle. 
However, as suggested by the results of this study 
one must expect that such field components are at 
least three times smaller than the field produced by 
flux sources. Thus, splitting this part from the total 
field is not trivial. It becomes clear in this context 
why there were always some doubts regarding the 
existence of a vortex field. 

The findings of this study show that vortex com¬ 
ponents will play a role in the genesis of MCG T- 
wave inversion, but they can not proof that closed 
loop currents yield the prominent contribution to the 
displacement of the MFM zero-line. From the 
simulated field patterns in Fig. 3. it can be seen, that 
the zero-line displacement occurs in the area close to 
the apex of the heart. The location of the apex 
changes from rest to stress due to the increased 
contraction of the heart. Thus, it can be suspected 
that also this spatial displacement of sources con¬ 
tributes to the MFM zero-line displacement during 
stress. It can also be observed that the isopotential¬ 
lines in the region closest to the apex are of concave 
shape at rest and of convex shape during stress. 
Therefore, the source displacement during stress is 
reflected also by the potential pattern. 

The location of the isopotential-zero-line showed 
only minor displacements in all simulations. This 
can be explained as follows: transmural hetero¬ 
geneity of repolarization leads to positive charges at 


the epicardium and to negative charges in the mid¬ 
wall (M-cells). Thus, the location of the electric 
potential zero-line is associated with the location of 
the heart base (here, the double layer of positive and 
negative charges is opened). 

The simulated zero-line displacement was of op¬ 
posite direction with respect to the measured. In this 
context it should be stated that the effect of stress on 
the source distribution, i.e. the fiber orientation and 
the twist on the heart could not be taken into account 
in this model. It is likely that such effects have a 
strong influence on the computed dipole moments. 
A MFM-zero-line displacement however, occurred 
only for the anisotropic model of the left ventricle. 
This supports the hypothesis that closed loop 
currents (which should not be mistaken with reentry 
currents) and a displacement of volumetric charges 
lead to MFM T-wave inversion. 
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